Two atoms in an anisotropic harmonic trap 
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We consider the system of two interacting atoms confined in axially symmetric harmonic trap. 
Within the pseudopotential approximation, we solve the Schrodinger equation exactly, discussing 
the limits of quasi-one and quasi-two-dimensional geometries. Finally, we discuss the application of 
an energy-dependent pseudopotential, which allows to extend the validity of our results to the case 
of tight traps and large scattering lengths. 
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PACS numbers: 34.50.-s, 32.80.Pj 

Atomic interactions at ultralow temperatures are of 
central importance for recent research on quantum de- 
generate gases 1] . A typical feature of experiments on 
ultracold matter is the presence of a weak trapping poten- 
tial, which modifies the properties of the cloud of atoms, 
while it does not affect the collisions of individual parti- 
cles. Development of optical lattice technology, however, 
has created systems where the atoms are tightly confined 
in the wells of optical potential Q{. In addition, the ex- 
perimental achievement of the Mott insulator phase Q 
has allowed for a precise control over a number of atoms 
stored in a single well. This has opened a way for exper- 
imental studies of interactions of individual atoms in the 
presence of trapping potential and, together with other 
approaches to micromanipulation of neutral atoms like 
atom chips 0, or tight dipole traps [|| , it represents 
a major candidate for the implementation of quantum 
information processing. A theoretical understanding of 
the dynamics of few atoms in deformed tight-confining 
geometries would be of great help in all these contexts. 

From the theoretical side, the analytical solution for 
two atoms interacting in a harmonic trap is known only 
for the spherically symmetric case 0,0. The correspond- 
ing problem for axially symmetric trap was studied nu- 
merically in |9(. However, there the authors considered 
only the limiting regimes of quasi-one and quasi-two- 
dimensional traps. In this letter we present the exact 
solution for the axially symmetric harmonic trap of ar- 
bitrary geometry. In particular, when the ratio of axial 
to radial trapping frequency is an integer, or the inverse 
of an integer, we give the explicit analytic form of the 
exact solution. In the other cases, we derive an efficient 
recurrence relation that allows for evaluating it. Further- 
more, we study the asymptotic behavior of eigenenergies 
and eigenfunctions in the limit of quasi-one and quasi- 
two-dimensional traps. 

A standard treatment of ultracold atom interactions is 
based on the replacement of a real physical potential by 
an s-wave delta-function pseudopotential. To extend the 
validity of this model interaction to the case of tight traps 
and large scattering lengths, one can utilize the concept 
of an effective, energy-dependent scattering length |10j . 
We discuss this idea and show how our results can be 
generalized to the case of magnetically tunable Feshbach 



resonances. 

We consider two interacting atoms of mass m confined 
in an axially symmetric harmonic trap with frequencies 
luj_ and lu z . In the following we use dimcnsionlcss vari- 
ables, in which all lengths are expressed in units of har- 
monic oscillator length a z = y/h/mu> z , and all energies 
are expressed in units of frw z . In these units the trap- 
ping potential is Vy(r) = ^(i] 2 p 2 + z 2 ), where t] = oj±/uj z 
and p 2 = x 2 + y 2 . We assume that the range of the in- 
teratomic potential is much smaller than the oscillator 
lengths a z and a± = y/h/mu>±, which guarentees that 
the interatomic potential is not distorted by the harmonic 
trap. For sufficiently low energies, the scattering is purely 
of s-wave type and we model the atom-atom interaction 
by a Fermi pseudopotential V(r) = 47ra5(r)J^r with s- 
wave scattering length a For the harmonic confining 
potential, the total Hamiltonian 



H = - + Vr(ri) + V T (r 2 ) + V(n - r 2 ), (I) 

can be splitted into center of mass part: i?cM = 
~hy 2 R + Vt(K), and the relative motion part: H le \ = 
-iV 2 . + V T (r) + V(V2r), where r = (n - r 2 )/v / 2 and 
R = (ri + r 2 )/v / 2- To solve the Schrodinger equation 
for the relative motion, we decompose the wave function 
in the basis of eigenstates of the noninteracting prob- 
lem, substitute this decomposition into the Schrodinger 
equation, and then extract the expansion coefficients by 
projecting onto noninteracting states Q. This yields the 
wave function of m z = states, with vanishing angular 
momentum along z-axis 
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The harmonic oscillator states with m z =/= vanish at 
r = 0, and they are not influenced by the pseudopo- 
tential. Eq. J2Jl represents the wave function which is 
not normalized, and is related with the single particle 
Green function of the anisotropic harmonic oscillator by 
'l'(r) = — 2G(r, 0). We note that the integral representa- 
tion © is valid for energies smaller than the ground state 
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energy of the harmonic oscillator: Eq = 1/2 + 77. The va- 
lidity of Eq. J2J), however, can be extended for E > Eq 
by means of the analytic continuation. 

The presence of the trapping potential implies the dis- 
crete character of the energy spectrum. The allowed 
values of energy E has to be determined from equa- 
tion: — 1/ (v2Va) = [{d/dr)r^(r)] r=0 , which results 
from derivation of Eq. (J2J, and expresses a boundary con- 
dition imposed by zero-range interaction. Investigation 
of the integral in Eq. |3J) for small values of r, shows 
that ^(r) behaves like l/(27rr) as r — > 0. This diver- 
gence is removed by the regularization operator (d/dr)r 
in the Fermi pseudopotential. Subtracting from the inte- 
gral J2J, the part which gives rise to the 1/r singularity, 
the condition for the eigenenergies can be rewritten as 



2*/a = F(-(E-Eo)/2,n), 



where 



T(x, rf) = / d/j 
Jo 



rje 



Vl - e"' (1 - e~"*) 



1 



(3) 



(4) 



For particular values of the anisotropy parameter 77, the 
function T(x, rf) can be calculated analytically. In the 
case of cigar shaped traps with 77 = n, where n is a posi- 
tive integer, we obtain 



!F(x, n) 



^T{x)J2 F(l,x;x+ \;e l ~* , 
m =i 2^/ttT(x) 

r(x + |) r(x-|)' 

(5) 

where F(a,b; c; x) denotes the hypergeometric function 
and r(x) is the Euler gamma function. It can be easily 
verified that the sum in Eq. JSJ involving complex roots 
of unity is a real number for x £ R. On the other hand, 
for pancake shaped traps with anisotropy parameter 77 = 
1/n, the following result holds 



T(x + m/n) 



n ^— ' r(x — 1/2 + m/n) 

m=0 y I'll 



(6) 



For n = 1, we recover obviously the well known re- 
sult for the spherically symmetric trap: T(x, 1) = 
-2-v/tt r(x)/r(^ - 1/2) 7]. We note that Eqs. © and 
© are derived from the integral representation (JIJ ap- 
plicable for x > 0, however, their validity for a; < is 
extended by virtue of the analytic continuation. 

In the general case, when 77 does not meet the con- 
ditions of Eqs. (J5J and the energy spectrum can be 
determined numerically. For E < Eq the function T(x, 77) 
is given by by Eq. while for E > Eq, one can utilize 
the following recurrence relation 

F(x, ??) - T(x + n, n) = r 1 ^T(x)/T(x + 1/2), (7) 

which can be easily derived from the definition of T(x, rf). 

From the practical point of view, the use of the exact 
results of Eqs. (J5J and © is efficient as long as n is not 



too large. To determine the energy levels in the limit 
of quasi-one- and quasi-two-dimensional traps, we derive 
the asymptotic form of T(x, rf) for 77 3> 1 and 77 <C 1. 

Let us first focus on the case of 77 ^> 1. Performing an 
expansion in the integral Q for large 77 and making use 
of the recurrence formula Q we arrive at 

F(x,rf) ^}{((±,l + x/ V ) + ^r(x)/r(x+±)}, (8) 

where Q(s, a) denotes the Hurwitz zeta function. This 
asymptotic formula is valid for x > —77, which corre- 
sponds to the range of energies E < Eq + 2rj. For the 
lowest excited states: < E — Eq <C 2?7 we approximate 
C(l/2, 1 + x/rf) by ({1/2,1) in Eq. ©, and match the 
resulting energy spectrum with the energy spectrum of 
two atoms in a one-dimensional trap. The latter is deter- 
mined by V2a 1D = T ((E - E)/2) /V ((E + 1 - E)/2) 
0- The two spectra are identical, provided that the 
one-dimensional scattering length is a±u = —1/rja — 
C(l/2, l)/\f2r\, which agrees with the value of the re- 
normalized scattering leng th derived for a quasi-one- 
dimensional waveguide |l3( . On other hand, for energies 



E < Eq, we can use T(x,rf) ~ y/wq ((1/2, x/rf), which 
follows from Eqs. (7J and (JSJ. This approximation, sub- 
stituted into J3J, leads to the condition determining the 
energy of a bound state: 



V2/a + (1/2, (Eq - £)/(277)) = 0, 



(9) 



which is identical to the known result derived for the 
quasi-one-dimensional waveguide 0, [w| . 

In the case of quasi-two-dimensional traps: 77 <C 1, we 
obtain the following approximate formula for T(x, rf): 



F(x,rf) 



77<1 



-<&(sc) - log(f7) - ip (x/rf) , (10) 



where 



$(x) = 2 - log(l + x) 
(2k)\ 



(11) 
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and ij}(z) = (d/dz)\ogT(z) denotes the digamma func- 
tion. This result is valid for x > — 1, which corre- 
sponds to energies E < E + 2. For the lowest ex- 
cited states: < E — Eq <C 2, we approximate Q(x) by 
$(0) in Eq. ifTUI) . and compare the resulting energy spec- 
trum to that of the two-dimensional system. In the two- 
dimensional trap, the eigenenergies of two interacting 
atoms are given by — log(2a2 D 7y) = -0 ((Eq — E)/(2rf)) 
|18| . In this way we find the value of the two-dimensional 
scattering length 0,20 for which both spectra are the 
same: a 2 D = exp[i(2? — \[2~k / 'azT>)}/\/2, where V = 
$(0) ~ 1.938. This result agrees with the value of 02D 
derived for a quasi-two-dimensional system without con- 
finement in the radial direction |15| . In the range of en- 
ergies corresponding to a bound state, we use an asymp- 
totic expansion of ip(x/rf) for x/r\ 3> 1 in Eq. i|ll)|) . which 
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be determined from the following expansions 
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FIG. 1: Energy spectrum of two atoms interacting via 
regularized delta potential in a three-dimensional trap with 
r\ — u)_i_/lu z = 100 (a) and r\ — 0.01 (b). Panel (a): The exact 
energy levels (solid lines) are compared with the energy spec- 
trum of the one-dimensional system with renormalized scat- 
tering length (dashed lines), and with the energies of a bound 
state calculated from Eq. @. (dotted lines - almost indistin- 
guishable from the solid ones). Panel (b): The exact energy 
levels (solid lines) are compared with the energy spectrum 
of the two-dimensional system with renormalized scattering 
length (dashed lines) , and with the energies of a bound state 
calculated from Eq. 1121 (dotted lines) . The scattering length 
a is scaled in h.o. units a z = yjhj- 



mui z - 



yields !F(x, rj) sw —$>(x)—logx. Substituting this approxi- 
mation into (j3J , we obtain the equation which determines 
the energy of a bound-state in quasi-two-dimensional 
traps: 



2ir/a = - E)/2) + log((£ - E)/2). (12) 

For a shallow bound-state {Eq — E <C 1) one can ap- 
proximate Q((E — E)/2) by $(0) and in this regime the 
binding energy is given by Eq — E — 0.288 exp(v / 27r/a) 

m 

Fig. n shows the energy spectrum of two interact- 
ing atoms calculated for r\ = 100 (a) and r\ = 0.01 
(b). Fig. IHa) compares the exact energy levels given 
by Eqs. J3J and JSJ, with the energy spectrum of 
the one-dimensional system with renormalized scatter- 
ing length aiD, and with bound-state energies calculated 
from Eq. Fig.^b) presents the exact result of Eqs. 
© and © , the energy spectrum of the two-dimensional 
system with renormalized scattering length a,2T>, and 
bound-state energy calculated from Eq. (|T2")) . We have 
not included the energy levels calculated from approxi- 
mations © and IjlOfl . which for 77 = 100 and 77 = 0.01 
are indistinguishable from the exact result. We observe 
that for E > Eq the one- and the two-dimensional spec- 
tra fit very well the exact eigenenergies, whereas they are 
incorrect with respect to the bound-state energies. The 
latter, however, are well described by Eqs. © and (|T2")) . 

We now turn to the calculation of wave functions. 
While for E < Eq they can be evaluated from the in- 
tegral representation |(SJ|, in the general case, they can 
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Here 8 = E — Eq, L m (x) and Hk(z) are respectively 
the Lagucrre and Hcrmite polynomials, D v (x) is the 
parabolic cylinder function and U(a, 6, z) denotes the 
confluent hypergeometric function. As it can be easily 
observed, the first expansion involves the harmonic os- 
cillator wave functions in the radial direction and the 
one-dimensional solution for two interacting atoms in the 
axial direction. We have verified that for elongated traps 
(?7 ^ 1), the first term of this series provides already a 
quite good approximation for the wave function of the 
lowest excited states. A similar feature is observed for 
the second series in Eq. I|14|) in the traps with 77 <C 1. 
Conversely, for energies E < Eq the two series involve 
generally several terms. In this regime, we can analyze 
the behavior of the wave functions on the basis of the 
integral representation © . Due to the complicated form 
of the latter integral, we focus here only on the limiting 
case of quasi-one- or quasi-two-dimensional traps, and in- 
vestigate only the behavior of the axial (p = 0) and the 
radial (z = 0) profiles of the wave functions. 

Expanding the integral in Eq. J2J) for 77 ^ 1, we obtain 
the axial (^> z (z) = *(zz)) and radial (^j_(p) = ^{pp)) 
profiles of the wave function, applicable for E < Eq 
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For |z|V— £ 3> 1, the main contribution to the sum in 
Eq. (|15fl comes from the first term. In this case the wave 
function exhibits the exponential decay, which is simi- 
lar to the behavior of one-dimensional bound-state in a 
free space: ^(z) ~ exp(— V— 2£|z|). On the other hand, 
the wave function in the radial direction has a Gaussian 
profile, characteristic for the ground-state of harmonic 
oscillator, whereas the divergent term l/(2irp) arises due 
to the interaction potential. 

In quasi- two-dimensional traps, for energies E < Eq, 
we found the following radial and axial profiles of the 
wave- functions 
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FIG. 2: The axial (p = 0) and the radial (z = 0) profiles of 
the ground-state wave function for two atoms interacting via 
a regularized delta potential with a = ±oo. The atoms are 
confined in a harmonic trap with rj = ui±/uj z = 100 (a), and 
r\ — 0.01 (b). The exact profiles (solid lines) are compared 
with the approximate results of Eqs. 115H - 117|I (dashed lines). 
All lengths are scaled to a z = \J h/mui z - 

where K (x) is a modified Bessel function. The asymp- 
totic behavior of Kq(x) for x 3> 1 is governed by Kq(x) ~ 
y/n/2xe~ x . Hence, for p\J — £ ^> 1, the sum in H7I) is 
dominated by the first term, and the asymptotic decay 
of the wave function in the radial direction is similar to 
the one observed for a bound state in two dimensions: 
^f(p) ~ Kq{\J—2£p). Along the tightly confined, ax- 
ial direction, the wave function has a Gaussian profile, 
which is modified at short distances by the interaction 
potential. 

The behavior of the ground-state wave function in the 
unitarity limit (a = ±oo) in the quasi-one-dimensional 
(ji = 100) and quasi-two-dimensional traps (rj = 0.01) 
is presented in Fig. [21 The figure compares the exact 
profiles evaluated from Eqs. (|13|l and H14() with the ap- 
proximate results of Eqs. (|15fl - l|17fl . We observe that all 
approximate curves fit quite well the exact functions. 



Finally we would like to stress that our derivation can 
be easily supplemented to include an energy-dependent 
scattering length [H [Hill. This extends the validity of 
the pseudopotential approximation to scattering lengths 
much larger than the trap size, and allows to properly 
describe the entire molecular spectrum. The energy- 
dependent effective scattering length is defined through 
the s-wave phase shift <5o: a s{E) = — t&n5o(k)/k, where 
hk is the relative momentum |l9j. The application of 
this model interaction in our derivations leads to substi- 
tution of a by a e g(E) in Eq. @ determining the eigenen- 
ergies, and requires a self-consistent solving for the value 
of E. For magnetically tunable Feshbach resonances, the 
s-wave phase shift is known analytically ^1 > an d in this 
case one can derive an explicit formula for a e g(E) |10| . 

In summary, we solved analytically the problem of two 
atoms interacting in an axially symmetric harmonic trap 
with arbitrary trap anisotropy. For integer ratios of the 
trapping frequencies we gave closed formulas for the so- 
lutions. Furthermore, by introducin g an effective energy 
dependence in the scattering length [Hj, [HI , we can find 
the solutions for any value of the latter. Therefore our 
result allows for a direct exact evaluation of the dynam- 
ics of a pair of interacting neutral atoms in very tight 
traps, possibly in reduced dimensionality and under an 
arbitrary external magnetic field, even in the presence 
of Feshbach resonances. Applications include a signifi- 
cant range of situations involving quantum control at the 
atomic level, from single-atom interferometry to quantum 
information processing. 
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